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Oscillators under the Minimal Length Uncertainty Relation 
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We analyze the position and momentum uncertainties of the energy eigenstates of the harmonic 
oscillator in the context of a deformed quantum mechanics, namely, that in which the commutator 
between the position and momentum operators is given by {x,p] — ih{l + Pfi^). This deformed 
commutation relation leads to the minimal length uncertainty relation Ax > (/i/2)(l/Ap + /3Ap), 
which implies that Aa; ~ 1/Ap at small Ap while As ~ Ap at large Ap. We find that the 
uncertainties of the energy eigenstates of the normal harmonic oscillator (m > 0), derived in Ref. [1], 
only populate the Aa; ~ 1/Ap branch. The other branch. Ax ~ Ap, is found to be populated 
by the energy eigenstates of the 'inverted' harmonic oscillator (m < 0). The Hilbert space in the 
'inverted' case admits an infinite ladder of positive energy eigenstates provided that Axmin = h^/p > 
\/2[t'-? /k\m\ ]^^'^. Correspondence with the classical limit is also discussed. 

PACS numbers: 03.65.-w,03.65.Ge 
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I. INTRODUCTION 

The consequences of deforming the canonical commu- 
tation relation between the position and momentum op- 
erators to 



(1) 



have been studied in various contexts by many authors 
[1-8]. The main motivation behind this was to use such 
deformed quantum mechanical systems as models which 
obey the minimal length uncertainty relation (MLUR) 
[9]: 



(2) 



This relation is expected on fairly generic grounds in 
quantum gravity [10, 11], and has been observed in per- 
turbative string theory [12]. The MLUR implies the ex- 
istence of a minimal length 



(3) 



below which the uncertainty in position, Aa;, cannot 
be reduced. For quantum gravity, Axmin would be the 
Planck scale, £p — \J HG n/<^, while for string theory this 
would be the string length scale, is = \fo! . The investi- 
gation of said model systems could be expected to shed 
some light on the nature of these, and other, theories 
which possess a minimal length scale. 

Note that the uncertainty in position which saturates 
the MLUR bound behaves as Aa; ^ 1/Ap for Ap < 
l/\fP, while Ax ~ Ap for l/\fP < Ap, as illustrated 




FIG. 1: The lower bound on Aa; as a function of Ap when 
they obey the minimal length uncertainty relation, Eq. (2), 
is shown in red. The blue line indicates Aa; — (fi/2)(l/Ap) 
while the green line indicates Aa; — {hp/2)Ap. 



in Fig. 1. While we are familiar with the Aa; ^ 1/Ap be- 
havior from canonical quantum mechanics, the Aa; ~ Ap 
behavior is quite novel. It behooves us to understand 
how it can come about, and how it can coexist with the 
canonical Aa; ^1/ Ap behavior within a single quantum 
mechanical system. To this end, we calculate Aa; and Ap 
for the energy eigenstates of the harmonic oscillator. 



2 2m 



(4) 
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the wave-functions of which were derived explicitly in 
Ref. [1]. We find that all of the eigenstates of this 
Hamiltonian inhabit the Aa; ~ 1/Ap branch, and not 
the other branch, as long as both the spring constant k 
and the mass m remain positive. The cross-over happens 
when the inverse of the mass, l/m, is allowed to decrease 
through zero into the negative, in which case all the en- 
ergy eigenstates will move smoothly over to the Aa; ~ Ap 
branch. The objective of this paper is to provide a de- 
tailed account of this result. 
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In the following, we solve the Schrodinger equation for 

the above Hamiltonian without assuming a specific sign 
for the mass m. The spring constant k is kept positive 
throughout. We find that the 'inverted' harmonic oscil- 
lator with fc > and m < admits an infinite ladder of 
normalizable positive energy eigenstates provided that 



where 



k\m\ 



1/4 



(5) 



(6) 



is the characteristic length scale of the harmonic oscilla- 
tor. The uncertainties Aa; and Ap are calculated for the 
energy eigenstates, and the above mentioned cross-over 
through 1/m — from the Ax ^ 1/Ap branch to the 
Aa; ~ Ap branch is demonstrated. 

We then take the classical limit of our deformed com- 
mutation relation and work out the evolution of the clas- 
sical harmonic oscillator for both the positive and nega- 
tive mass cases. It is found that for the 'inverted' m < 
case, the time it takes for the particle to travel from 
X = — oo to a; = +oo is finite, demanding the compactifi- 
cation of a;-space, and also rendering the classical prob- 
ability of finding the particle near the origin finite. This 
provides a classical explanation of why 'bound' states 
are possible for the 'inverted' harmonic oscillator in this 
modified mechanics. 



II. QUANTUM STATES AND UNCERTAINTIES 

A. Eigenvalues and Eigenstates 

The position and momentum operators obeying Eq. (1) 
can be represented in momentum space by [13] 

p = p. 

The inner product between two states is 

dp 



(7) 



(1 + /?p2) 



This definition ensures the symmetricity of the operator 
X. The Schrodinger equation for the harmonic oscillator 
in this representation is thus 



*(p) = E^{p) . (9) 



Here, we do not assume m > as usual, so the kinetic 
energy term can contribute with either sign. A change of 
variable from p to 



P = arctan(y^p) 



(10) 



maps the region — oo < p < oo to 

and casts the x and p operators into the forms 

dp 

P = -^tan(^p) , 
with inner product given by 



(11) 



(12) 



(13) 



Note that x is the wave-number operator in p-space, so 
the Fourier coefficients of the wave-function in p-space 
will provide the probability amplitudes for a discretized 
x-space. Eq. (9) is thus transformed into: 



h^k 1 ^ r- 



*(p) = E^{p) . (14) 



In /9-space, the potential energy term kx'^/2 effectively 
becomes the kinetic energy, and the kinetic energy term 
p^/2m effectively becomes a tangent-squared potential 
which is 'inverted' when 1/m < 0. We next introduce 
dimensionless parameters and a dimensionless variable 
by 



2\m\Ej3 

c = ^ 

s — ) 



Aa;i, 



(15) 



where the length-scale a was introduced in Eq. (6). The 
dimensionless variable ^ is in the range 



TT TT 

~2k ^ 2k 



(8) and the inner product is 



_K_ r/'" 

a/^ J-w/2k 



(16) 



(17) 



The dimension of the inner product has all been absorbed 
into the prefactor l/\/^- The Schrodinger equation be- 
comes 



1 



■ tan + e 



*(0 = , (18) 



where the minus sign in front of the tangent-squared po- 
tential is for the case m > 0, and the plus sign for the 
case m < 0. Let ^(^) = f{s), where s = sinK^, 
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FIG. 2: Plots of A+ (blue) and A_ (red), Eq. (23), as functions 

of = Phy^k]^. 



c = cosK^ = Vl — s^, and A is a constant to be de- 
termined. The variable s is in the range 



-1< s < 1 
with inner product given by 

The equation for /(s) is 

(l_s2)/"-(2A + l)s/' 



0*(s)V(s) 



(19) 



(20) 



+ 



/ = o 



(21) 



We fix A by requiring the coefficient of the tangent 
squared term to vanish: 



A(A-1)T A 



= 



The solutions are 



A = 



(m > 0) 
(to < 0) 



(22) 



(23) 



where we have chosen the branches for which A > 1/2 
to prevent the inner-product, Eq. (20), from blowing up 
at the domain boundaries. Note that 1 < A4. while i < 
1 in the limit 



A- < 1, and that A± 



l3hy/k\^\-> 



00. The dependence of A± on k is shown in Fig. 2. The 
A_ branch does not extend below = 2. Setting A = A+ 
for TO > 0, and A = A_ for m < simplifies Eq. (21) to 

(1-s2)/"-(2A + 1)s/'+(^-A)/ = 0, (24) 

the sign of m being encoded in the value of A. Since /(s) 
should be non-singular at s = ±1, we demand a polyno- 
mial solution to Eq. (24). This requirement imposes the 



following condition on the coefficient of /: 



(25) 



where n is a non-negative integer [14]. Eq. (24) becomes 

(1 - s^)f" -{2X + l)sf' + n{n + 2\)f = 0, (26) 

the solution of which is given by the Gegenbauer polyno- 
mial: 



(27) 



The Gegenbauer polynomials satisfy the following or- 
thogonality relation: 



27rr(n-f 2A) 
[2^r(A)]2n! (n + X) 



Sn 



(28) 

The energy eigenvalues follow from the condition 
Eq. (25). Replacing A with A±, we find 



r(±) - 



[n^ + {2n+ l)Ad 
+ {2n + 1)A± 



A± A± 



K^[n^+n+-] + (2n+ 1) 



± 1 



(29) 



or in the original dimensionful units. 



1 



-f (2n-h 1)A 



± 



2/3|m 



A±|A±-1 

1\ /^2^2^2 



± 1 



1\ (3\m\haj' 
^^2) 2 



-) v/(Aa:„,in)4 ± 4a4 



(30) 



where uj — ^kl\m\. For the to > O case, we can take the 
limit Axmin = ~^ Oi ^-i^d we recover 



lim = kc? 



= hw 



(31) 



For the to < O case, it is clear that we must have Axmin ^ 
^/2a for the square-root in Eq. (30) to remain real. This 
is the condition we cited in Eq. (5). Therefore, the limit 
Axmin = ^a/? — > cauuot be taken in this case for non- 
zero a. The two cases converge when \m\ — oo, at which 
a — 0, and we find that the energy levels in that limit 
are 

, lim = |(A.T„i„)2(n + l)2. (32) 

|m|— >-cx3 Z 
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Thus, the 1/m > and 1/m < cases connect smoothly 
at 1/m — 0. 



by: 



where 




n! (n + A) 



27rr(n + 2A) 
1 



(34) 



The wave-functions for the first few energy eigenstates for 
several representative values of A are shown in Figs. 3. 



B. Expectation Values and Uncertainties 

Using the wave-functions derived above, and the for- 
mula provided in the appendix, the expectation values of 
X, p, i^, and for the energy eigenstates are found to 
be 



(n, A| X |n, A) — , 
(n, A| p |n, A) — , 



(2A-1) 



{n,X\x^\n,X) = {h^l3)- 



i,X\p \n,X) = - 



1 f2n + l 



P V2A- 1, 
giving the uncertainties in (a;) and (p) as 



(35) 



(A + n)[(2A- l)?i + A 



(2A-1) 



1 / 2n+l 



(36) 



For fixed n and fixed /3, Ax„ is a monotonically decreas- 
ing function of A in the range ^ < A < 1, and a mono- 
tonically increasing one in the range 1 < A. Ap„ on the 
other hand is a monotonically decreasing function of A 
throughout. Eliminating A from the above expressions, 
we find 



1 



> 



2V?AP« 
1 / 1 



2 VV?Ap„ 



1 + (3Apl {2n + 1)2 + (3Apl 
+ ^Ap^) , (37) 



where the equality in the second line is saturated for the 
n = case only. The first line gives the curve on the 
Ap-Ax plane that the point {Apn, Axn) follows as A is 
varied. Differentiating with respect to Apn we find 

Axn 



d{Apn) 



W/3 



thus given 






\ A=100 


(33) 
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FIG. 3: A-dependence of the wave-functions of the first few 
energy eigenstates. The A > 1 values correspond to m > 0, 
while the | < A < 1 values correspond to m < 0. The A = 1 
case corresponds to the limit |m| — >■ oo. 



P^Api~{2n+lf 



1 + pApl 



(2n + 1)2 + pApl 



(38) 



indicating that the curve is flat at the A = 1 point where 
Apn — \J {2n + l)//3 and Axn reaches its minimum of 
Axi-ain{n + 1). Therefore, the A = 1 point is the turn- 
around point where the uncertainties switch from the 
Aa; ~ 1/Ap behavior to the Aa; ~ Ap behavior. To 
go from one branch to another one must flip the sign of 
the mass m. 

Eliminating n from Eq. (36), we find 



Aa;n 

AXini 



(l + ^Apl) 1 + (2A - l)2/3Ap2 . (39) 
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FIG. 4: Ap versus Aa; for the lowest six energy eigenstates 
of the harmonic osciUator. Aa; is in units of Axmin = fty^i 
while Ap is in units of ~ ft/Axmin- The location of the 

state along the curves shown is determined by the value of A 
defined in Eq. (23). The A = 1 points correspond to the case 
1/m = 0. As 1/m is increased to the positive side, the value of 
A will increase away from one and the state will move toward 
the left along the Aa; ~ 1/Ap branch. If 1/m is decreased into 
the negative, the value of A will decrease toward 1/2, and the 
state will move toward the right along the Ax ~ Ap branch. 
The n — curve (shown in red) corresponds to the minimal 
length uncertainty bound, Eq. (2). 



This gives the curve on which the points (Ap„, Aa;„) fall 
on for constant A. In particular, for A = 1 this reduces 
to 



1 + /?Ap2 



(40) 



and gives the 1/m = boundary between the 1/m > 
and 1/m < regions in Ax-Ap space. These properties 
of the uncertainties have been plotted in Fig. 4. 



1. A — ^ cx) 

The (3 limit only exist for the m > case where 
A = A+. As /3 — >■ 0, the parameter diverges to infinity 
as 



A = A. 



->■ oo , 



(41) 



where w = yjk/m. In that limit, the Gegenbauer poly- 
nomials become Hermite polynomials: 



lim n\\-'^'^C^{x/^/\) = Hn{x) . 

A— >oo 

Noting that as /3 0, we have 
s ~ 



we can conclude that 



lim n!A-"/2(^;t(s) 

A— >-oo 



V Xmfujj 



Similarly, 



lim = lim ( 1 — . ^ ^ — ) 



A— >c 



A— >c 



,2 \ 

Xmhio J 



exp 



(42) 



(43) 



(44) 



2mhuj 



(45) 
(46) 



Using Stirling's formula 

the normalization constant can be shown to converge to 

1 1 



\n/2 

lini ^ 

A-s-oo ni 



•\/2"n! \/mhjj'K 



(47) 



Therefore, 



lim ^>\^\p) 



1 



1 



•\/2"n! mhtoTT 



exp 



\ jj^ ( P \ 
2mhw) "vVmfio;/ 



(48) 



which are just the usual harmonic oscillator wave- 
functions in momentum space. The energy eigenvalues 
reduce to the usual ones given in Eq. (31). The uncer- 
tainties reduce to the usual ones as well 



C. Limiting Cases /s^x 



As shown in Fig. 4, the value of A determines where 
the uncertainties (Ap„, Ax„) are along their trajectories 
given by Eq. (37), with A > 1 keeping the uncertainties 
on the Aa; ~ 1/Ap branch of the trajectory, while \ < which satisfy 
A < 1 keeping them on the Aa; ^ Ap branch. Let us 
consider a few limiting values of A to see the behavior of 
the solutions there. 



Ap„ 



X^oo 



1 



/ hmuj 



AxnApn ^^°°> h{n+^j > I 



(49) 



(50) 
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2. \^\ 

The limit A = 1 is reached when /3 and k are kept 
constant while |m| is taken to infinity. When A = 1, the 
Gcgcnbaucr polynomials become the Chebycheff polyno- 
mials of the second kind: 



where 



sin(n +1)9 



UnicOsO) = 



smf 



while the normalization constant reduces to 



Therefore, 



-cU„{s) . 



(51) 



(52) 



(53) 



(54) 



The orthonormality relation for the Chebycheff polyno- 
mials is 



j y/l-X^ Unix) Um (x) dx = 



(55) 



and we can see that the correct normalization constant is 

obtained. Since the argument in our case is s = sin y^p, 
it is more convenient to express the Chebycheff polyno- 
mials as 



U^nisinO) = {-!) 



^cos[{2n + 1)0] 

cos 6 



U2n+i{sme) = (-1) 



^ ( ^.„ sin[(2n + 2)g] 



cos^ 



(56) 



for n 



0, 1, 2, • • • . This will allow us to write 



2n+l 



i-iy 



cos 



(2n + l)s/pp 



(57) 



(-l)"\/-sin (2n + 2)V?p 



The energy eigenvalues in this limit were given in 
Eq. (32). Here, our procedure of keeping the spring con- 
stant k fixed while taking |m| to infinity maintains the 
finitcncss of £"„ , while taking the kinetic energy contribu- 
tion to En to zero. From the n-dependcnce of the ener- 
gies, we can see that, in this limit, the problem reduces to 
that of an infinite square well potential, of width tt/v^, 
in /9-space. Indeed, the effective potential in p-space was 







DO 



for p ^ ± 
at p = ± 



2y?' 



(58) 



This can also be seen from the form of the energy eigen- 

functions, which have reduced to simple sines and cosines. 
We will see in the next section that the classical solution 
also behaves as that of a particle in an infinite square 
well potential in p-space in the same limit. 
The uncertainties become 



AXn 



A->-l 



Aa;niin('^ + l) 



. A^i {2n+l) 



(59) 



as was shown in Fig. 4. Note that 



1 



Aprr 



n+1 



l3Apn = Ax J,,,,, 

y'zn + 1 

< AXmin{n + l) = AXn , 



(60) 



the bound being saturated only for the ground state n = 
0. 



A^i 



The A — > i limit is reached as Axy^in — >■ V^a when 
m < 0. In this limit, the Gegenbauer polynomials be- 
come the Legendre polynomials, Cn^^{s) = Pn{s), while 
the normalization constant reduces to 



The wavefunctions are 

^'h^^\p) ^ /2n + l 



^/iPnis) 



(61) 



(62) 



Note that the orthonormality relation for the Legendre 
polynomials is 



Pn{x) Pni{x) dx = + l ^""^ ' 



so these wave-functions are properly normalized. The in- 
tegrals for (i;^) and (p^) diverge in this limit, so both Ax„ 
and Apn are divergent for all n. However, the energy, 
which is the difference between k{x^)/2 and (p^)/2|m|, 
stays finite: 



E- 



III. 



CLASSICAL STATES AND 
UNCERTAINTIES 



(64) 



As we have seen above, for values of j3 which maintain 
the inequality Axmin = > \/2a, the harmonic os- 

cillator Hamiltonian admits an infinite ladder of positive 
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energy eigenstates even when m < 0. Furthermore, these 
are states with finite Aa; and Ap, implying that the par- 
ticle is 'bound' close to the phase space origin, just as in 
the m > case. But how can a particle be 'bound' for an 
'inverted' harmonic oscillator? To gain insight into this 
question, we solve the corresponding classical equation of 
motion. 



A. The Classical Equations of Motion 

We assume that the classical limit of our commutation 
relation, Eq. (1), is obtained by the usual correspondence 
between commutators and Poisson brackets: 



in 

Therefore, we assume 



B] ^ {A,B} 



{x,x} = , 
{P,P} - 0, 
{x,p} - + 



(65) 



(66) 



Then, the equations of motion for the harmonic oscillator 
with Hamiltonian given by 



1 
2m 



P 



(67) 



X = {x, H} = -{l + (3p^)p, 
m 

p = {p,H} = -k{l + Pp^)x. (68) 

We allow m to take on either sign: if m > 0, then x 
and p will have the same sign; if m < they will have 
opposite sign. Note that, even though the equations of 
motion of x and p have changed, the total energy will 
still be conserved. Consequently, the time-evolution of x 
and p in phase space will be along the trajectory given 
hy H — constant. For the m > case this will be an 
ellipse, while for the m < case this will be a hyperbola. 

To solve these equations, we change the variable p to p, 
which was introduced in Eq. (10) for the quantum case. 
Then, the equations become 



1 



—kx . 



tan(v/?p) 



cos2(^p) 



1 



2m/3 dp 



ia.v?{^p) 



(69) 



Therefore, 



-kx 



2m(3 dp 



tan2(V^p)l , (70) 



which integrates to 



k 

m/3 



tan2(y^p) -C 



(71) 



where C is the integration constant. Since we must have 
p^ > 0, the range of allowed values of C will depend on 
whether to > or m < 0. We will consider the two cases 
separately. 
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FIG. 5: The dependence of the classical behavior of a pos- 
itive mass particle in a harmonic oscillator potential on the 
parameter C = 2mi5/3 = /3m^cj^y4^, where E is the particle's 
energy, and A is the amplitude of the oscillation in x. The 
undeformed /3 = case is shown in red. The other four cases 
are C — 1 (orange), C = 2 (green), C = 4 (blue), and C — 8 
(purple). 



B. m > case 

When TO > 0, we introduce the angular frequency 

T 



as usual. Then, Eq. (71) becomes 



/3p' = uj' 



C-tan2(y^p) 



and taking the square-root, we obtain 



(72) 



(73) 



(74) 



In this case, we must have C > for the content of 
the square-root to be positive. Separating variables, we 
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FIG. 6: The classical behavior of a positive mass particle in a 
harmonic oscillator potential with modified Poisson brackets, 
Eq. (66), in the limit m — >• oo with C = 2mEl3 where E 
and P are kept fixed, pit) and x{t) take on the behavior of 
the position and momentum of a particle in an infinite square 
well. 



obtain 



The left-hand side integrates to 

__Sdp__ 

VC-tan2(V?p)_ 
1 



— ztiodt 



(75) 



Therefore, 

= arcsin 



C 



l + C 



sin(v/^p) 



sin |±\/1 + Cuj{t - to) I 



(76) 



(77) 



where to is the integration constant. Without loss of gen- 
erality, we can choose the sign inside the curly brackets 
to be plus. Setting the clock so that to = 0, we obtain: 

\/ fSmk x(t) = —^^-^p 



y/C{l + C) cos(VTTCwt) 
y/l + C'cos2(^/rTCwt) 
^p{t) = tan(y^p) 



^Jl + C cos^ {■k/TTC Lot) ' 
and the energy is given by 



(78) 



\/ I3mk x{t) 
> 0. 



2/3m 
C 

2/3m 



2/3m 



(79) 



The period of oscillation T is no longer equal to 2tt/lu 
when /3 ^ 0. It is now 



Let 



Then 



T = 



A 



271 1 



/ C 

j3m,k 



E = -kA^ 
2 



(80) 



(81) 



(82) 



and we can identify A as the oscillation amplitude in x. 
If we take the limit (3-^0 while keeping A fixed, we find: 



x{t) 
Pit) 



—Acos{ujt) , 
Amuj sin(ajt) 



(83) 



which shows that the canonical behavior is recovered in 
this limit. The behavior of the solution when /? ^ is 
compared with the /3 = limit for several representative 
values of C in Fig. 5. 

Another interesting limit is obtained by setting C = 
2'mEP and letting m oo while keeping E and /3 fixed. 
In that limit, 



and we find 



% ^/2E(3k = ^/^kA = ujo , (84) 



x{t)/A 



sin(ajoi) 



m^Qo sin(a;oi) 
)■ ' 



I cos(a;ot)| 
cos(ajot) 

I cos(woOI 



(85) 



The behavior of the solution in this limit is shown in 
Fig. 6. The motion of a particle in an infinite square well 
potential (in p-space) is reproduced, in correspondence 
to the quantum A — > 1 limit. 
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C. m < Case 



For the to < case, by an abuse of notation, let us set 






2 

7T 










nr 
V \m\ ■ 


(86) 






4 








Then, Eq. (71) becomes 




3 


2 


/ / n 








= ^2 


tan^ {y/l3p) - C 


(87) 






/ /~ 4 
1 ^ 









The integration constant C can have either sign in this 
case. We will consider the three cases C < 0, C > 0, and 
C — separately. 



1. C < (positive energy) case 
For the C < case, the square-root of Eq. (87) gives us 



y/Pp = ±ujJtan^{^p) + \C\ 



VP dp 



Therefore, 



tan2(V^p) + Id 

The left-hand side integrates to 

VPdp 
Vtan2(V^p) + \C\ 
1 



± u) dt 



(89) 



sm{y/Pp) 

1 



sin(V^Py 



Therefore, 
sin(v^p(i)) 



: arcsinh 



sm 



> 1) 

(1^^1 = 1) 
{\c\ < 1) 

(90) 
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FIG. 7: The dependence of the classical behavior of a neg- 
ative mass particle in a harmonic oscillator potential on the 
parameter C — -~2pEl3 — — j3\m\kB'^ , where E is the parti- 
cle's energy, and B is the distance of closest approach to the 
origin in a;-space. Note that due to the negative mass, p{t) 
is negative when x(t) is positive, and vice versa. The unde- 
formed /3 = case is shown in red. The other four cases are 
C — —-^ (orange), C = — -j (green), C = — 1 (blue), and 
C = -4 (purple). 



\c\ 



l-\C\ 



sinh 



±yT^[q^(<-to)J {\c\<i) 

(91) 



where to is the integration constant, which we will set to 
zero in the following. From this, we find: 

a/ /3|TO|fc x(t) ~ — 



^ Vic|(|q-i) cos{±V\c\^ut) ^1^1 ^ 



|C|cos2(±^|C| - lujt) - 1 
1 



{\c\^i) 



^ V|C|(1 - \C\) cosh(±yT^c.t) ^1^1^^^ 
'l - |C|cosh2(±y"l -\C\ujt) 
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\C\ sin(±^|C| -lujt) 



± 



UJ t 



\C\ sinh(±^l- |C|a;t) 
_ Jl-\C\cosh^{±y^l-\C\ujt) 



(|C|>1) 

i\c\^i) 

(|C| < 1) 
(92) 



In all three cases, we have 



2 ^ ' 2\m\ 
1 



2p\m\ 

\C\ 
2l3\m\ 



^Jl3\m\k x{t) 
> 0. 



2/3|m| 



Let 



B = 



/ \c\ 

f3\m\k 



Then 



E = -kB^ 
2 



(94) 



(95) 



and we can identify B as the distance of closest approach 
to the origin (aka impact parameter). Taking the hmit 
/3 — >■ while keeping B fixed, we find: 



x(t) 
Pit) 



=Fi? cosh(±wt) , 
B\m\uj sin]i{±Lut) , 



(96) 



which recovers the canonical solution. This behavior of 
x{t) and p{t) for the /? = case is compared with that in 
the (3 case for several representative values of C in 
Fig. 7. 

It should be noted that for any finite value of C < 0, 
it only takes a finite amount of time for the particle to 
get from {x,p) = (±oo,±oo) to {x,p) — (±oo,^cx)), or 
equivalently, for y/^p to evolve from =F7''/2 to ±tt/2. We 
will call this time T/2 for reasons that will become clear 
later. T is given by: 



1 



1 



T ^ - X < 



^/\C\^ 
1 
1 



arccosh 



This dependence on C < is shown in Fig. 9 



(|C| > 1) 

(iq = l) 
{\C\ < 1) 



(97) 
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FIG. 8: The dependence of the classical behavior of a neg- 
ative mass particle in a harmonic oscillator potential on the 
parameter C — —2\m\Ej3 — /Sp^jjn, where E is the particle's 
energy, and pmin is the momentum of the particle at the origin 
a; = 0. The undeformed P = Q case is shown in red. The other 
four cases are C — + (orange), C = +\ (green), C = +1 
(blue), and C — +4 (purple). 



2. C > (negative energy) case 



For the C > case, taking the square-root of Eq. (87) 
yields 



(98) 
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Therefore, 



^tan2(V^p) -C 
The left-hand side integrates to 

Vtan2(V^p)-C 
1 



±U! dt 



(99) 



VT+C 



arccosh 



1 + C 
C ' 



(100) 



Therefore, 

sin(V^p(0) 

= ± 



C 



1 + C 



cosh 



±VT+c uj{t - to 



(101) 



where to is the integration constant, which we will set 
to zero in the following. The sign on the argument of 
the hyperbolic cosine is also irrelevant so we will set it to 
plus. From this, we find: 

VP . 



Vl3\m\kx{t) = - — 



^C(l + C) smh{y/TTCLut) 



1 - Csmh^{^/TTCuJt) 

tan(V£p) 

C cosh(\/rTCa;t) 



± 



y^l - Csmh^{VTTCujt) 



(102) 



and 



E 



■x{tf 



2/3|to| 



Pit? 
2\m\ 

Vl3\m\kx{t) 



2l3\m\ 




FIG. 9: LoT as a function of —C — 2|m|£'/3, where oj — 
^Jk/\m\, and E is the particle's energy. T/2 is the time it 
takes for the particle to traverse the entire trajectory. 



C 

'2(3\rr 



Let 



Then 



< 



Pmir 



E 



'2\m\ ' 



(103) 



(104) 



(105) 



and we can identify Pmin as the magnitude of the momen- 
tum that the particle has at the origin x = 0. Taking the 
limit /3 — > while keeping Pmin constant, we find 



x[t) > =F 1 — — smh(±u;t) 



Pit) 



> ± Pmin cosh(±a;t) . 



(106) 



As in the C < 0, the time T/2 it takes for the particle to 
travel from x — ^oo to x — ±oo is finite. T is given by: 



4 .1 

T = . arcsmh — ;= 

LOy/TTC Vc 



(107) 



This dependence on C > is also shown in Fig. 9. 

3. C = (zero energy) case 

For C = 0, Eq. (87) leads to 

^/Pp = ± cj tan(y^p) , 



or 



y/pdp 



— zL ujdt , 



tan {V0p) 
which can be integrated easily to yield 



In 



sm 



(VPp) - ±uj{t~to) 



or 



sin(V^p) = ± e 



±tj(t-to) 



(108) 



(109) 



(110) 



(111) 



with all combinations of signs allowed. Set the clock so 
that to = 0. The solutions for the t > region are 



a/ (3\m\k x{t) = —^^p 



±- 



1 



y/Pp{t) = tan{^/^p) = ± 



(112) 



Ve2"* - 1 



The particle starts out at {x,p) = (±cx), ±cx)) and asymp- 
totically approaches the origin, taking an infinite amount 
of time to get there. This behavior is show in Fig. 10. 
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FIG. 10: The classical behavior of a zero-energy particle with 
negative mass in a harmonic oscillator potential. The particle 
starts out at a; = — oo at time t = and asymptotically 
approaches the origin. 



FIG. 11: The periodic behavior of the negative-mass parti- 
cle in a harmonic oscillator potential in compactified a;-space. 
The example shown is for C — —1. 



D. Compactification 

As we have seen, when m < and /3 7^ 0, it only takes 
a finite amount of time for the particle to traverse the 
entire classical trajectory as long as the energy of the 
particle is non-zero. This means that we must specify 
what happens to the particle after it reaches infinity. For 
this, we could either compactify x-space so that the par- 
ticle which reaches x — ±00 will return from x = ^00, 
in which case the momentum of the particle will bounce 
back from infinite walls at p = ±00, or we could compact- 
ify p-space so that the particle which reaches p — ±00 will 
return from p = ^00, in which case the position of the 
particle will bounce back from infinite walls at x — ±00. 

Here, we choose to compactify x-space so that the 
|m| — >■ c» Hmit of the m < 0, C < solution will match 
the m 00 hmit of the m > 0, C > solution. This 
choice also agrees with the boundary condition we im- 



posed in the quantum case, in which the wave-function 
in p was demanded to vanish at the domain boundaries 
p = ±7r/2-v^, which corresponds to placing infinite po- 
tential walls there. The m > 0, C > solution was given 
by Eq. (85). Taking the \m\ ^ 00 limit of Eq. (92) while 
keeping B fixed, we find 



and 



(113) 



y^p(i) 
x{t)/B 



|m|- 



arcsm 



|m|- 



|m|- 



sin(±cjoi) 

I COs{±ljjQt)\ 

cos(±woi) 
I cos(±a;oi)| 



(114) 



which formally agrees with Ec^. (85), and if graphed will 
lead to a figure similar to Fig. 6. The one significant 
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difference is, however, that when the particle jumps from 
X — ±A to X = =fA in the m > case it goes through 
x = 0, while when it jumps from x = ±B to x = =p_B in 
the m < case, it must go through a; = oo. 

By compactifying x-space, all motion when E will 
become oscillatory through a; = oo, and the T calculated 
above becomes the oscillatory period. As an example, we 
plot the x-compactified solution for C = — 1 in Fig. 11, 
for which the period is T = A/uj. Note that the period 
for the TO < 0, C < solution in the limit of |m| — > oo 
becomes 



lim T 

|m|— >oo 



lim 



1 



the arccosine providing a 7r/2. 



: arccos 



27r 

(115) 



E. Classical Probablities 

Consider the to < 0, C = —2\m\EI3 < case. y^P 
evolves from ~tt/2 to tt/2 in time T/2, that is: 

-7r/2v? P 

(116) 



T 
2 



T/4 

dt = 

-T/4 



— dp 



--n^ ^tan2(V^p) + 
Thus, we can identify 

2 V? 
'^^ /tan2(V?p) + |C| 



(117) 



as the classical probability density of the particle in p- 
space. The classical probablity in p- and x-spaces can be 
defined in a similar manner: 



Pip) = P{p) 



dp 
dp 



P{x) 



ojT (1+ I3p2)^/\C\TW ' 

T 
2 



^P\m\kx^ - |C|(/3|TO|fc2a;2 _ \c\ 



|C|(a;2-S2) + S2 



(118) 



These probability functions are plotted in Fig. 12 for the 
case C = -2\m\EP = -1. 

Comparing the energies of the quantum and classical 
solutions given in Eqs. (30) and (93), we can conclude 
that the correspondence is given by the relation 



-C = 



+ (2n + 1)A 
A(l-A) 



i<A<l 



(119) 
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FIG. 12: Classical probabilities in p-, p-, and a;-spaces of a 
negative mass particle in a harmonic oscillator potential for 
the case C = —2\m\EI3 — —1. The time dependence of this 
solution was shown in Fig. 11. Though the trajectory of the 
particle is not confined to a finite region of phase space, the 
classical probabilities of finding the particle near the phase- 
space-origin is still finite due to the finiteness of the time it 
takes for the particle to traverse the entire trajectory. 



We expect the quantum and classical probabilities to 
match for large n. As an example, we take A = | and 
n = 30, which correspond to: 



C = 



B 



-5044 

^^min 



^A(l-A) ^ 



1.52, 



30Ax„ 



(120) 



The comparison of the quantum and classical probabil- 
ities for this case in p-, p-, and x-spaces are shown in 
Fig. 13. If we average out the bumps in the quantum 
case, it is clear that the distributions agree, up to the typ- 
ical quantum mechanical phenomenon of seepage of the 
probability into energetically forbidden regions. Thus, 
the existence of 'bound' states with a finite Aa; and Ap 
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FIG. 13: Comparison of quantum and classical probabilities 
in p-, p-, and a;-spaces for a negative mass particle in a har- 
monic oscillator potential for the case A = I and n = 30, 
which corresponds to the classical C — — 2|m|iJ/3 = —5044. 
The quantum distribution in a;-space is discrete due to the 
existence of the minimal length. There is also some seepage 
of the probability into classically forbidden regions in a;-space 
as is expected of quantum probabilities. 



in the quantum case can be associated with the fact that 
the particle spends a finite amount of time near the phase 
space origin in the classical limit. 



IV. SUMMARY AND DISCUSSION 

We have solved for the eigenstates of the harmonic 
oscillator hamiltonian under the assumption of the de- 
formed commutation relation between x and p as given 
in Eq. (1), with the objective of calculating their uncer- 
tainties in position and momentum. 

For the normal harmonic oscillator with positive mass 
{1/m > 0), the eigenstates are found on the Ax ^ 1/Ap 
branch of the MLUR, where decreasing 1/m leads to 
larger Ap, and thus smaller Ax. Somewhat surpris- 
ingly, 1/m can be decreased through zero into the neg- 
ative, thereby 'inverting' the harmonic oscillator, while 
still maintaining an infinite ladder of positive energy 



eigenstates so long as the condition Axmin/v^ > a = 
[ ?i^/fc|m| ]^/^ is satisfied. There, the eigenstates are 
found on the Ace ^ Ap branch of the MLUR, where 
decreasing 1/m away from zero further into the negative 
leads to larger Ap, and thus larger Aa;, with both diverg- 
ing as a approaches the above bound from below. The 
1/771 = line separating the Ax ^ l/Ap and Aa; ^ Ap 
regions is given by Eq. (40). 

Taking the classical limit by replacing our deformed 
commutator with a deformed Poisson bracket, we solve 
the corresponding classical equations of motion and find 
that the solutions for the 'inverted' harmonic oscillator 
are such that the particle only takes a finite amount of 
time to traverse its entire trajectory. This leads to a fi- 
nite classical probability density of finding the particle 
near the phase space origin, and provides and explana- 
tion of why 'bound' states with discrete energy levels are 
possible in the quantum case. 

One significant difference between the classical and 
quantum cases is, however, that the classical system has 
no restriction on the sign of the energy, whereas the quan- 
tum system only allows for positive energy eigenstates. 
The latter is guaranteed by the above mentioned condi- 
tion on Axniin and a. Indeed, the condition is equivalent 
to 



E = ^(Axf 



1 



2|m 



■{Apf > 



(121) 



when one assumes Ax ^ {h(3/2)Ap. 

We have found that the energy eigenstates of the har- 
monic oscillator only populate the Ax ~ 1/Ap branch 
of the MLUR when the mass is positive, and only the 
Aa; ^ Ap branch when the mass is negative. A natu- 
ral question to ask is whether some superposition of the 
energy eigenstates could cross over the 1/m = line to 
the other side. For the free particle, which can be con- 
sidered the k — limit of the harmonic oscillator, the 
answer is in the affirmative. In that case, the uncertain- 
ties of the energy eigenstates are Ap = and Ax = 00. 
Superpositions of these states, with uncertainties posi- 
tioned anywhere along the MLUR bound, can be easily 
constructed, as will be discussed in detail in a subsequent 
paper [15]. Is a similar 'cross-over' possible for the fc ^ 
case? 

Another interesting question is whether it is possible 
to construct classically behaving 'coherent states' in ei- 
ther of these mass sectors, and whether their uncertain- 
ties are contained as they evolve in time. In particu- 
lar, when the mass is negative, the classical equations 
of motion calls for the particle to move at arbitrarily 
large speeds. What is the corresponding quantum phe- 
nomenon? We are also assuming that Eq. (1) embodies 
the non-relativistic limit of some relativistic theory with 
a minimal length. From that perspective, the infinite 
speed that the negative-mass particle attains seems prob- 
lematic. Would such states not exist in the relativistic 
theory, or will they metamorphose into imaginary mass 
tachyons? These, and various other related questions will 
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be addressed in future works. 



we can write 
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Appendix A: An Integral Formula for Gegenbauer 
Polynomials 

The following integral formula is necessary in calculat- 
ing the expectation value of . For two non- negative 
integers m and n such that ni < n, and A > i, we have 



= < 



47rr(m + 2A) 
(2A-l)[2^r(A)]^m! 



I 



if n — m = even, 

if n — m = odd, 

(Al) 



where c = Vl — s^. We were unable to find this result 

in any of the standard tables of integrals [14] though 
Mathematica seems to be aware of it. Here we present a 
proof. 

We start from recursion relations which can be found 
in Ref. [14]. In section 8.933 of Gradshteyn and Ryzhik, 
we have: 



2A 



c^+\x)-xc::ilix) 



nC^ix) = 2x\xC^tl{x)-C^:tlix) 



(A2) 



Eliminating C^^l{x) and then shifting A by one unit, we 
obtain 



A- 1 



A- 1 



Cr^(x) , (A3) 



which is Equation 22.7.23 of Abramowitz and Stegun. 
Iterating this relation, we deduce that 



k 



A - 1 



Ck+ii^) = C^{x)+j2(^^)c^^~+\i^)■ (A4) 



Since 



C^ix) = 1 ^ Ct'ix) 

C^{x) = 2Xx = (j^ ] Ct\x) , (A5) 



i=0 ^ ^ 

cL+lix) = E (f^r) • (A6) 

i=0 ^ ' 

Thus, the even C^'s can be expressed as a sum of the even 
C^"^'s, and the odd C^'s as a sum of the odd C^~^'s. 
Invoking the orthogonality relation, Eq. (28), which is 
valid when A > — |, it is clear that 



/: 



(A7) 



for A > |. So for the integral of Eq. (Al) to be non- 
zero, m and n must be both even, or both odd. For two 
non-negative integers k and i such that k < i,we find 



„2A-3 



E ( ^^fV ) "^'^ 

. i=0 

E 



ds 



k I 

EE 

i=0 3=0 



2i + A - 1\ /2j + A - 1 
A- 1 J \ A- 1 

"\^^-^C^-\s)C^-\s)ds 



1 



^ /2i + A- 1^ ^ 



i=0 



A- 1 



27rr(2i + 2A-2) 



(2i + A- l)[2^-ir(A- 1)] (2i)! 



2tt 



[2^-'mYt. 



(2i + A-l)r(2i + 2A-2) 



(2i)! 



„2A-3 



k t 



i=0 j=0 ^ ' ^ 



ds 



_ ^ /2i + A^ ^ 



j=0 



X f^c''-'C^r^\{s)C^f+\{s)ds 
2'KT{2i + 2X-l) 



A- 1 
27r 



(2i + A)[2^-ir(A-l)]'(2i + l)! 



[2^-^r{x)]-' to 



{2i + X)T{2i + 2X-l) 
{2iTT)\ 



(A8) 
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The sums in the above expressions are given by 



In particular, the diagonal elements are given by 



E 

i=0 
k 

E 

i=0 



(2i + A- l)r(2i + 2A-2) 
(2i + A) r(2i + 2A - 1) 



{2i + l)\ 



T{2k + 2A) 
~ 2(2A- l)(2fc)! ' 

r(2fc + 2A + 1) 
2(2A- l)(2fc + l)! ■ 
(A9) 



These relations can be proved by induction in k. Putting 
everything together, we obtain Eq. (Al). 

Using this formula, we find the matrix elements of the 
operator to be: 

{m\p^\n) = (n|p^|m) 



1 



+ 



2V(A + m)(A + n) / n!r(2A + m) 

2A- 1 y m!r(2A + n) 
for n — m = even, m < n 



for n — m = odd. 

(AlO) 



{n\p \n) = — 



1 /2n+ 1 



/? V2A- 1 



The expectation value of Sp' is obtained from 



(All) 



(A12) 
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